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|  This  report  presents  an  advanced  constitutive  model  to  describe  the  complex  behavior  ol 
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I  crophysically  bus.  i  constitutive  relationships  to  model  deformation  and  damage  processes  in  a 
j  ceramic.  The  total  strain  is  decomposed  into  clastic,  plastic,  and  microcracking  components 
:  The  model  parameters  for  AD85  ceramic  were  determined  using  the  data  from  spin  Hopkinson 
j  bar  (SHB)  and  bar-on-bar  experiments  under  uniaxial  stress  stale  and  plate  impact  experiment 
i  under  uniaxial  strain  state.  To  further  validate  the  generality  of  the  model  parameters,  model 
|  *ng  a  diagnostic  ballistic  experiment  in  which  a  steel  projectile  impacted  an  ADS  5  ceramic 
front-faced  thick  aluminum  plate  was  considered.  In  this  experiment,  stress  histories  were  mca 
sured  in  the  target  by  embedded  manganin  and  carbon  stress  gauges.  The  results  from  the  nu 
merical  simulations  of  the  ballistic  experiment  using  a  shock  wave  propagation  based  finite 
j  clement  code  successfully  matched  the  measured  stress  historv. 
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Introduction 


Advanced  armor  design  calculations  to  understand  and  evaluate  the  performance 
of  armor  elements  often  employ  hydrocodes  (shock  wave  propagation  based  finite 
element/dtfference  numerical  codes).  For  realistic  design  calculations  using  hydro¬ 
codes,  the  impact  behavior  of  armor  materials  must  he  accurately  described  with 
physically  based  constitutive/failure  models.  This  is  especially  true  for  advanced 
armor  materials  such  as  ceramics.  Since  constitutive  and  damage  models  lor  ceramics 
are  rarely  available  in  hydrocodes,  several  preliminary  versions  of  new  models  have 
been  implemented  into  different  hydrocodes.  Most  of  these  new  models  are  in  the 
process  of  being  validated  and  evaluated. 

The  high  velocity  impact  behavior  of  ceramics  is  dominated  by  stiffness  loss.  The 
inelastic  deformations  are  due  to  microcracking  and  dislocation  generated  nicroplastit 
flows.  In  the  microcracking  approach,  the  (strain)  response  of  a  single  crack  to 
an  externally  applied  stress  deld  is  calculated  using  the  appropriate  stress  free  hound 
ary  conditions  on  the  crack  surfaces.  The  derived  strain  components  include 
the  effects  of  microcracks  in  a  brittle  material  through  a  crack  density  parameter 
However,  the  bulk  material  itself  is  treated  as  a  continuum  so  that  all  the  stress 
strain  equations  can  be  derived  under  the  continuum  mechanics  based  theoretical 
framework.  The  microplasticity  in  a  brittle  solid  is  often  attributed  to  the  dislocation 
motions  in  the  vicinity  of  microflaw  tip  regions.  In  brittle  solids,  large  scale  grain 
distortions  are  usually  absent. 

Unlike  metals,  most  ceramics  pulverize  upon  high  velocity  impact  thereby  eliminating 
the  possibility  of  any  post-impact  microscopic  measurement/analysis.  The  American 
Physical  Society  conference  proceedings  [1]  report  several  microstructural  studies 
on  ceramic  targets  recovered  undet  low  velocity  (below  Hugoniot  Elastic  Limit!  im¬ 
pact  conditions.  The  microscopic  evaluations  of  the  recovered  targets  provide 
useful  information  on  microcrack  initiation  under  compressive  loading  conditions 
Since  the  effects  of  nucleation  and  growth  of  microcracks  on  the  constitutive  behaviors 
of  ceramics  are  important  in  the  modeling,  the  low  velocity  impact  studies  are  olten 
important. 

In  the  absence  of  any  post-impact  measurements  on  most  ceramic  targets, 
the  ceramic  behavior  under  high  velocity  (above  Hugoniot  Elastic  Limit)  impact  is 
usually  evaluated  through  stress  and  velocity  measurements  at  the  back  of  the  target 
plates.  These  measurements  yield  dynamic  properties  such  as  the  Hugoniot  stress  and 
spall  strength.  In  these  experiments,  the  ceramic  targets  are  backed  by  plastic  materi¬ 
als  such  as  plexiglass  and  lithium  fluoride.  Rajendran  and  Cook  |2|  conducted  an  in 
depth  literature  review  of  impact  behavior  of  ceramics  and  reported  the  need  tor 
ceramic  constitutive  and  damage  models  in  advanced  hydrocodes.  Subsequently, 
during  1988  through  1992  a  few  new  models  have  been  developed  and  reported  to 
describe  the  impact  behavior  of  ceramic  materials. 

Johnson  and  Holmquist  [3|  considered  a  phenomenological  approach  in  their 
model  formulation.  The  strength  variation  with  respect  to  pressure,  strain  rate,  and 
damage  was  modeled  using  a  set  of  empirical  relationships.  In  their  model,  the 
strains  are  due  to  elastic  and  plastic  deformations.  No  attempt  was  made  to 
characterize  the  inelastic  deformation  due  to  microcrack  opening  and/or  sliding. 


Damage  is  defined  as  the  ratio  of  accumulated  strain  and  a  fracture  strain.  This  deli 
nition  of  damage  is  similar  to  the  metal  fracture  model  of  Johnson  and  Coo  k  |4 1. 
Johnson  and  Holmquist  included  data  from  both  static  and  dynamic  tests  under  uniax¬ 
ial  compressive  loading  in  their  model  parameters  reduction  scheme. 

Steinberg  [5]  adopted  his  metal  model  equations  for  high  pressure  dynamic  load¬ 
ing  to  describe  the  impact  behavior  of  ceramics.  In  this  model,  the  compressive  yield 
strength  and  shear  modulus  are  varied  with  respect  to  pressure  and  temperature.  In 
addition,  the  yield  strength  is  multiplied  by  the  normalized  shear  modulus.  The  strain 
rate  effects  are  introduced  through  a  simple  power  law.  There  are  two  constants  in 
the  power  law.  The  two  constants  are  determined  using  the  fracture  toughness  depen 
dent  relationship  of  Grady  [6|  and  the  quasi-static  and  split  Hopkinson  bar  (SHB) 
data.  However,  the  tensile  damage  is  modeled  through  the  Cochran  and  Banner  spall 
model  [7],  Steinberg’s  model  assumes  no  compressive  damage  in  the  ceramics. 
Steinberg  successfully  reproduced  the  measured  velocity  histories  from  plate  impact 
experiments  for  various  ceramic  materials. 

Addessio  and  Johnson  [8]  presented  a  microphysical  model  to  describe  the  com¬ 
plex  behavior  of  ceramics.  In  their  model,  the  inelastic  strain  is  assumed  to  be  due 
to  microcracking  of  the  ceramics.  They  modeled  crack  opening  under  tensile  pressure 
loading  and  crack  sliding  under  compressive  pressure.  The  plastic  How  or  pore  col¬ 
lapse  effects  are  not  built  into  the  model;  however,  an  extension  of  this  model  to 
include  such  effects  is  feasible.  Damage  is  described  by  a  crack  density  parameter. 
The  damage  evolution  is  described  through  energy-balance  based  failure  surfaces. 
Addessio  and  Johnson  introduced  arbitrary  modifications  to  the  model  and  successfully 
reproduced  the  measured  velocity  histories  for  several  ceramics  under  one-dimensional 
strain  conditions.  This  three-dimensional,  continuum  mechanics  model  has  not  been 
validated  under  any  other  stress-strain  states. 

Rajendran  and  Kroupa  [9|  presented  a  ceramic  model  in  which  they  assumed  an 
elastic-viscoplastic  behavior  for  compressive  loading  and  microcracking  behavior  for 
tensile  loading.  In  general,  attributing  the  inelastic  strains  entirely  to  dislocation- 
based  plastic  flow  is  not  appropriate  for  brittle  ceramics.  However,  this  model  may 
be  useful  io  describe  crushing/pulverization  of  the  ceramic  materials  under  high  pres¬ 
sures  (>10  GPa)  and  dominantly  tensile  damage.  The  tensile  cracking  is  described 
using  the  rock  fragmentation  relationships  derived  by  Grady  and  Kipp  |1()|.  Since  it 
is  a  fragmentation  based  tensile  fracture  model,  the  model’s  applicability  is  limited 
and  any  extension  of  this  model  includes  the  less  appealing  plastic  work  or  plastic 
strain-based  compressive  damage  description. 

Recently,  Rajendran  [11|  reported  a  continuum  mechanics  based,  three- 
dimensional  constitutive  model  to  describe  the  complex  behaviors  of  ceramic  materials. 
This  model  [12  through  1 4 J  is  based  on  microcrack  nucleation  and  growth,  as  well 
as  pore  collapse  mechanisms.  Damage  is  defined  in  terms  of  an  average  crack  den¬ 
sity  and  is  treated  as  an  internal  state  variable.  To  keep  the  model  formulation  rela¬ 
tively  simple,  the  damage  nucleation  is  not  modeled  and  the  microcracks  arc  assumed 
to  be  present  prior  to  loading.  This  scalar  damage  model  incorporates  the  effects  of 
different  damage  processes  under  tension  and  compression  using  fracture  mechanics 
based  fracture  criteria. 


f 


The  stiffness  reduction  due  to  microcracking  is  modeled  using  Margolin’s  [15] 
analytically  derived  damaged  moduli.  The  damage  evolution  under  both  tensile  and 
compressive  loadings  is  formulated  based  on  a  generalized  Griffith  criterion.  The 
pore  collapse  effects  are  modeled  using  viscoplastic  equations  derived  from  Gurson's 
pressure  dependent  yield  function  ( 1 6 ] .  This  model  decomposes  total  strains  into 
elastic  and  dislocation  motion-controlled  plastic  components.  The  elastic  strains  are 
further  decomposed  into  average  strains  in  the  intact  matrix  material  and  strains 
due  to  crack  opening/sliding.  The  strain  rate  and  pressure  effects  on  strength  arc 
implicitly  described  by  the  degraded  elastic  moduli  based  stress-strain  equations 
Rajendran  and  Grove  [14]  and  Grove,  et  ai.  [12]  successfully  applied  the  model  to 
describe  the  impact  behavior  of  AD85  under  uniaxial  stress  and  one-dimensional 
strain  conditions.  This  three-dimensional,  continuum  damage  model  ]11|  has  been 
incorporated  into  the  EPIC-2  finite  element  code  ]17],  Model  constants  were  deter 
mined  for  the  AD85  ceramic  using  static,  SHB.  plate  impact,  and  bar-on-bar  impact 
test  data.  This  report  describes  the  application  of  the  Rajendran-Grove  ceramic- 
model  to  a  problem  in  which  a  steel  projectile  impacts  a  layered  ceramic  target.  The 
main  objectives  are:  (1)  to  establish  generality  of  the  model  constants,  and  (2)  to 
demonstrate  the  model’s  ability  to  predict  the  measured  stress  history  under  muitiaxial 
loading  conditions.  The  center  portion  of  the  target  in  a  ballistic  test  initially  experi¬ 
ences  one-dimensional  strain,  and  later  muitiaxial  strains  due  to  the  release  waves.  • 

Only  the  shock  and  incipient  damage  phases  of  the  penetration  process  are  modeled. 

Vincent  and  Chang  [18]  performed  a  ballistic  experiment  in  which  a  steel  projectile 
impacted  a  layered  AD85  (aluminum  oxide)  target.  In  the  experiment,  the  shock 
stresses  were  measured  using  embedded  manganin  stress  gauges  at  two  locations  inside 
the  target.  The  measured  stress  history  is  influenced  by  the  various  shock  release  ^ 

waves:  elastic-plastic  compressive  shock  waves,  release  waves  from  the  edges  o|  the 
projectile-target,  and  fracture  waves  due  »o  ceramic  fracturing.  Therefore  the  model 
ing  of  a  ballistic  experiment  under  such  complex  wave  interactions  is  indeed  useful 
in  the  ceramic  model  validation.  In  the  "Rajendran-Grove  Ceramic  Model"  sec¬ 
tion.  the  salient  features  of  the  ceramic  modei  are  briefly  described.  The  model  pa¬ 
rameters  determination  scheme  for  AD85  ceramic  is  also  outlined.  The  "Modeling  * 

the  Ballistic  Impact  Pressure  Measurements"  section  briefly  describes  the  ballistic  exper 
iment  and  discusses  the  various  features  of  the  ceramic  model  in  reproducing  the  ex¬ 
perimentally  measured  stress  history.  The  "Summary  and  Conclusions"  section 
follows  the  above  sections.  The  damaged  moduli  expressions  of  Margolin  (15)  are 
provided  in  Appendix  A. 


Rajendran-Grove  Ceramic  Model 

Conventionally,  the  impact  behavior  of  a  material  is  described  through  a  strength 
model  and  an  cquation-of-state  (EOS).  The  strength  model  describes  the  variation  of 
strength  with  respect  to  the  strain  rate,  temperature,  and  pressure.  The  strength  is 
expressed  through  the  von  Mises  stress  vTJt  ,  where  Ji  is  the  second  invariant  of 
the  stress  deviators.  Therefore,  the  strength  mode1  involves  the  calculation  of  the 
deviatoric  part  (S}J)  of  the  stress  tensor  The  bulk  (volumetric)  behavior  of  the 
material  is  described  by  the  EOS.  This  involves  the  calculation  of  the  pressure 
(mean  stress)  part  of  the  stress  tensor. 
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Constitutive  Relationships 

The  total  strain  is  decomposed  into  elastic  and  plastic  strains  as 


e  a  =  e  ..  +  e 


where  the  elastic  strain  consists  of  the  elastic  strain  of  the  intact  matrix  material 
(e™),  and  the  strain  due  to  crack  opening/sliding  ( f  ,c, ) . 


e  _  m  .  c 

f  ij  f  >j  »j  ' 


Since  the  strain  components  due  to  microcracking  are  elastic  ( e  ^ ) .  represents  both 
crack  opening  and  closing.  The  matrix  elastic  strain  and  the  microcracking  strains  arc 
both  proportional  to  the  applied  stress  field.  The  plastic  strains  are  calculated  from 
viscoplastic  flow  equations.  In  the  Rajendran-Grove  ceramic  model,  Johnson  and 
Cook  [4)  type  strength  model  is  employed  to  describe  the  strain  rate  dependent 
strength  behavior.  When  voids  (pores)  are  present  in  the  ceramic,  the  strain  compo¬ 
nents  due  to  pore  collapse  are  calculated  from  Gurson’s  pressure-dependent  plastic 
flow  equations  (16).  The  elastic  strains  in  the  microcracked  void-free  aggregate  mate¬ 
rial  are  related  to  the  stress  tensor  as  follows: 


£  a  “  Qjkl  ^kl  • 


where  C,^  are  the  effective  compliance  tensor  components  of  the  microcracked 
material.  If  Cjjio  is  analytically  inverted  to  the  stiffness  tensor  components  Ml|kl.  the 
resulting  stress  state  is 


aij  ~  M,jk|  £|^j 


The  elements  of  the  stiffness  matrix  C,,^.  derived  by  Margolin  [15.19]  and  Budianskv 
and  O’Connell  [20),  are  described  in  Appendix  A.  1  he  total  stress  is  decom¬ 
posed  into  deviatoric  stress  and  pressure  components: 


a,j  ~  Su  +  P4, 


The  elastic  stress-strain  relationship  between  the  deviatoric  stresses  and  the  cone 
sponding  deviatoric  strains  is  given  by 


2  G  (e, 


4 


where  e,,  are  the  total  deviatoric  strains,  e  P  are  the  devialonc  plastic  strains,  and  G 
is  the  degraded  shear  modulus  due  to  microcracking  and  pores  When  there  are  no 
pores.  G  will  correspond  to  G  (see  Appendix  A)  which  is  the  shear  modulus  ol 
microcracked  ceramics.  Therctore,  the  effective  shear  modulus  (including  the  effects 
of  both  microcracks  and  pores)  G(=  R„G)  is  calculated  using  the  modified 
Mackenzie's  relationship  [2 1 .22 1: 


Ru  =  ( I  -  1") 


(6K.  +  12G  )f  j 
(9K  +  SG  )  j 


(  " ) 


» 


Rg  is  the  shear  modulus  correction  factor  tor  spherical  pore  content  in  the  ceramic 
f  is  the  porosity  content  (void  volume  fraction),  and  K  is  the  microcracks  degraded 
bulk  modulus.  The  equation  of  state  is  described  bv  the  following  modified  Mie- 
Gruneisen  relationship. 


where 


Ph  =  (K.K  +  ihu~  +  /T  x  // 1 )  i‘/i 

_  _V 

u  (  1  vo  *  is  the  elastic  volume  compressibilitv  (strain)  and  ,}>.  ,{■>.  and  ,it  are  the 

empirical  parameters,  E  is  the  Mie-Gruneisen  parameter.  /».,  is  ihc  matertal's  initial 
density.  1„  is  the  initial  value  of  internal  energy,  and  1  is  the  current  internal  cncrg\ 
The  Mackenzie's  correction  factor  R^  is  given  by 


Rk  = 


(1  ~ 


+ 


3K  r 

4G 


(  Mil 


» 


where,  Rk  K  is  the  degraded  bulk  modulus  due  to  pores  and  microcracks.  In  the 

absence  of  voids  and  microcracks.  Equation  S  reduces  to  the  Mie-Gruneisen  EOS  lot  » 

the  undamaged,  (lawless  material. 

Definition  of  Damage 

In  the  ceramic  model,  microcrack  damage  is  measured  in  terms  ol  a  dimensionless 
microcrack  density  y,  where  • 


» 


N  l 
‘  o  ‘  max 


N*  is  the  average  number  of  microflaws  per  unit  volume  in  the  ceramic  which  is  a 
model  constant.  amax  is  the  maximum  microcrack  sue  which  is  treated  is 
an  internal  state  variable.  Microcrack  extends  when  the  stress  state  satisfies  the 
Griffith  criterion  [23 J.  As  the  microcracks  extend,  y  increases  and  the  stresses  relax 
In  the  model  N*  is  assumed  to  be  constant  Therefore,  the  increase  in  y  is  due  to 
increase  in  the  crack  size  In  the  "Rajendran-Grove  Ceramic  Model"  and  Model 
ing  the  Ballistic  Impact  .  ressure  Measurement"  sections,  the  microcrack  growth  law 
is  defined  using  a  dynamic  fracture  mechanics  based  relationship. 

In  general,  the  phases  of  damage  evolution  consists  of:  (I)  nucleation  of 
microcracks  with  some  initial  crack  density,  (2)  growth  of  microcracks,  and  (3)  coales 
cenr  of  the  microcracks  at  some  critical  crack  size.  At  present,  no  means  exist  to 
experimentally  measure  and  quantify  the  variation  of  crack  size,  crack  density,  and 
number  of  cracks  with  respect  to  time.  However,  the  effects  of  microcracks  evolution 
can  be  introduced  into  the  constitutive  equations  through  appropriate  phvsicallx  based 
laws/equations. 

Microcrack  Nucleation 

One  simplified  approach  is  the  following:  The  microcracks  are  assumed  to  exist 
prior  to  any  loading:  therefore,  the  initial  crack  density  y0  is  determined  from  the  rela 
tionship  (see  Equation  11)  using  the  values  for  the  two  material  constants:  N*  and 


The  microdamage  (microcrack  size)  will  not  increase  until  the  generalized  Griffith 
criterion  (Margolin  |24|  and  Dienes  [ 25 1 )  is  satisfied  either  under  shear  loading  or 
under  tensile  loading.  Note  that  the  shear  loading  could  occur  under  both  tensile 
and  compressive  pressures.  The  microdamage  rate  is  zero  until  the  applied  strain 
energy  release  rate  G|  exceeds  a  critical  value  Gc  The  strain  energy  release  rate 
(note  that  the  rate  refers  to  crack  extension,  not  the  "time")  under  mode  1  crack 
opening  for  a  general  loading  (o:!  *  (I)  is  given  by: 


g|  = 


2  Ki-  + 


akk  + 


.  i^j^k 


where  <7kk  are  the  normal  stresses,  <7,k  and  a jk  arc  the  shear  stresses.  am;iv  is  the 
size  of  the  largest  microcrack,  and  v  and  E  are  undamaged  Poisson's  ratio  and 
Young's  modulus,  respectively.  The  repeated  index  ’k’  does  not  mean  summation. 
When  the  crack  surfaces  are  perpendicular  to  tensile  principal  stresses  o\,  then 
Equation  12  is  evaluated  with  akk  =  C7|  and  a Jk  =  a]k  =  1).  Under  general  loading 
conditions  (©^*0),  the  strain  energy  release  rate  for  mode  II  and  mode  III  is  ex 
pressed  as 
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where  r0  is  the  cohesion  stress  and  fx  is  the  friction  coefficient.  In  this  case,  only 
modes  II  and  III  are  active  and  the  normal  stresses  serve  only  to  resist  the  shearing 
stresses.  Defining  Gmax  to  be  the  maximum  of  all  values  of  Gj  and  G  ,  ,  micro¬ 
cracks  are  assumed  to  have  extended  if  Gmax  exceeds  the  critical  crack  energy  release 
rate  Gc.  where 


Gc  =  2  T 


and 


T  = 


Keff(1 
2  E 


(  15) 


Equation  15  defines  the  surface  tension  T  as  a  function  ol  effective  critical  frac¬ 
ture  toughness  (Kujf).  Poisson’s  ration  (r).  and  Youngs  modulus  (E)  in  the  undam¬ 
aged  material.  In  the  model  application,  the  critical  effective  stress  intensity  factor  is 
assumed  to  he  equal  to  simply  the  static  fracture  toughness  (K|k  )  under  crack  open¬ 
ing  mode.  At  present,  the  fracture  toughness  is  treated  as  a  material  constant  which 
will  not  vary  with  strain  rate.  John  [26]  experimentally  determined  the  strain  rate  de 
pendent  fracture  toughness  values  for  brittle  concrete.  Though  such  a  study  does  not 
exist  for  ceramic  materials,  especially  at  high  strain  rales,  it  is  still  possible  to  include 
strain  rate  effects  on  the  fracture  toughness  in  the  model  formulation  based  on  his 
study. 

Damage  Growth  Model 

The  damage  evolution  equation  is  derived  from  the  fracture  mechanics  based  rela¬ 
tionship  (27|  for  a  single  crack  propagation  under  dynamic  loading  conditions.  An 
evolution  law  for  the  state  variable  amax  is  described  through  a  strain  energy  release 
based  microcrack  growth  rate  law  of  the  form 
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(  16) 


where  Cr  is  the  Rayleigh  wave  speed.  Gc  is  the  critical  strain  energy  release  for 
microcrack  g/owth.  and  G\  is  the  applied  strain  energy  release.  The  model  constants 
nj  and  ni  can  be  used  to  limit  the  microcrack  growth  rate,  nj  and  n^  are  damage 
growth  constants.  Since  the  crack  growth  based  damage  rates  are  different  under 
tensile  and  compressive  loadings,  these  constants  will  be  assigned  different  values. 
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When  the  crack  density  reaches  a  critical  value  of  0.75  (see  Reference  i  1  ).  the 
model  assumes  that  the  microcracks  have  coalesced,  leading  to  pulvet l/ation  oi  the 
ceramic.  Henceforth,  the  material  has  no  strength  in  tension  and  its  compressive 
strength  follows  a  Mohr-Coulomb  law.  as  in 


Y  =  I  "D  +'3PP'PS>0  »hcn V  =  11.75 


where  Y  is  strength.  P  is  pressure,  and  ap  and  /3p  arc  model  constants  for  the  pulver¬ 
ized  material.  The  pressure  is  computed  from: 


ec  >  0 


-Knce  <  0 

l  V  v.  V 

where  ec  is  engineering  elastic  volumetric  strain  (see  Equation  2t  tind  Kp  is  the  hulk 
modulus  for  pulverized  material.  With  this  approach,  each  pulverized  element  in  a 
finite  element  mesh  may  have  its  own  distinct  values  of  Kp  and  Gp  .  It  is  also 
possible  to  input  new  values  for  K  and  G  as  bulk  and  shear  moduli  of  the  pulverized 
ceramic;  this  will  increase  the  number  of  mode!  constants  by  two. 

Model  Parameters  Determination 

The  static  compressive  strength  of  ADX5  ceramic  has  been  measured  at  about 
2  GPa.  The  SHB  test  data  (strain  rate  =10  s)  indicate  that  its  unconfined  dynamic 
compressive  strength  is  about  3  GPa  flij.  From  plate  impact  test  measurements,  the 
dynamic  compressive  strength  of  AD85  at  extremely  high  strain  rates  (  =  10^  s)  is 
about  5.5  GPa.  This  strength  increase  is  due  to  a  combined  strain  rale  and  confined 
pressure  effect.  While  the  confining  pressure  effect  was  implicitly  modeled  by  the 
stress-strain  relationship  (see  Equation  4).  the  compressive  strength  variation  with 
respect  to  strain  rate  is  described  bv  (Johnson  and  Cook  model  with  Cy  =  0): 


Y  =  C]  ( I  +  Cy  In  f  )  . 


The  constants  C]  and  are  calibrated  using  the  strength  values  obtained  from  the 
plate  impact  and  SHB  tests.  The  corresponding  values  are:  C\  =  4.0  GPa  and 
C-\  =  0.03  GPa.  When  the  strain  rates  tire  low'  (<1000/scc)  and  the  stress  state  is 
uniaxial,  the  ceramic  strength  is  controlled  by  the  microcracking  feature  of  the  model, 
therefore,  the  SHB  data  tire  reproduced  by  adjusting  the  microcrack  model  constants, 
especially  the  frictional  coefficient  u.  The  model  assumes  plastic  flows  in  ceramics 
when  the  shock  pressures  are  above  the  Hugoniot  clastic  limit.  The  total  strain 
includes  elastic,  cracking,  and  plastic  components.  The  plastic  strains  tire  calcu¬ 
lated  using  conventional  plastic  How  theories. 
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Thete  arc  six  constants  in  the  Rajendran-Grove  ceramic  model  to  describe  the 
microcracking  behavior:  Kk  ,,u.  N  * ,  ao .  ni.and  n^.  Note  that  the  cohesion  stress 
in  Equation  13  is  assumed  to  be  zero  and.  therefore,  this  constant  is  not  considered 
in  the  present  study.  The  preliminary  set  of  the  model  constants  is  determined  Iron; 
SHB  and  plate  impact  experimental  data. 

To  assure  generality  of  the  model  constants,  manganin  gauge  measured  stress  data 
from  the  bar-on-bar  impact  experimental  configuration  was  also  considered  The  best 
suitable  values  for  the  constants  are  obtained  based  on  the  ability  to  reproduce  the 
experimental  data  from  SHB,  plate  impact,  and  har-on-har  impact  configurations. 
Rajendran  [11]  and  Rajendran  and  Grove  [14|  reported  this  model  determination 
scheme  and  determined  the  model  constants  for  ADK5  ceramic  The  corresponding 
constants  are  given  in  Table  1. 


Table  1  Model  constants  for  AD85  ceramic 


Symbol 

Value 

Description 

Kic 

3  MPa  vm- 

Static  fracture  toughness 

ft 

0.72 

Coefficient  of  friction 

K 

1  83  x  10,0m'3 

Microcrack  density  (numbers.'volume) 

ao 

58  x  10  6  m 

Initial  microcrack  size 

n  2 

007 

Crack  growth  rate  power  index 

1  0 

Tensile  crack  growth  rate  index 

nT 

0  1 

Compressive  crack  growth  rate  index 

The  value  for  K|c  is  obtained  from  fracture  mechanics  handbooks.  With  lack  ol 
any  microscopic  measurements  a  large  value  for  the  number  of  flaws  in  the  ceramics 
N*  has  been  arbitrarily  assumed.  It  may  also  be  assumed  that  the  number  ol  llaws  is 
proportional  to  the  number  of  grains  per  unit  volume  and  attempt  to  estimate  this 
constant.  However,  in  the  present  work  nr)  such  attempt  was  made. 

The  plate  impact  simulations  revealed  that  the  initial  maximum  crack  size  a„  and 
the  crack  growth  indices  nj  and  n?  controlled  the  shape  of  ihe  unloading  portion  of 
the  stress-time  history  profile  of  the  plate  impact  data.  The  arrival  of  a  very  weak 
spall  signal  at  point  B  in  Figure  I  is  found  to  be  sensitive  to  the  crack  size.  In  the 
simulation,  the  shape  between  points  B  and  C  was  influenced  by  the  growth  con¬ 
stants.  Therefore,  the  initial  estimates  for  these  two  damage  growth  constants,  n  t 
and  n2,  were  calibrated  to  reproduce  the  stress  profile  measurements  in  a  bar 
(uniaxial  stress  state)  and  in  a  plate  (one-dimensional  strain  state). 

In  metals,  the  experimentally  measured  crack  propagation  speed  under  mode  1 
(crack  opening)  is  equal  to  the  fraction  of  the  Rayleigh  wave  speed  |27).  However, 
nj  is  set  ;  one  for  ADK5  ceramic.  The  crack  propagation  speed  under  shear  (mode  11 
or  mode  III)  is  relatively  lower  than  under  mode  I  and.  therefore,  a  value  of  0.1  is 
assumed  arbitrarily  for  nj.  For  nj,  a  value  of  0.07  was  determined  based  on  the 
model’s  ability  to  reproduce  the  rod-on-rod  experimental  data  (see  Figure  2). 


Figure  1 .  A  comparison  between  model  and  measured  stress  history 
in  a  plate  impact  experiment  on  AD85  target. 


The  model  parameters  estimation  scheme  requires  a  trial  and  error  basis 
of  adjusting  the  crack  growth  indices  between  different  experimental  configurations 
Unfortunately,  there  is  no  one  set  of  unique  values  for  the  impact  damage  model 
parameters.  This  is  true  for  all  the  ceramic  models  in  References  3,  5,  8.  9.  and  11. 
At  best,  a  suitable  set  of  values  can  be  successfully  determined  through  a  trial  and 
error  basis  of  reproducing  a  variety  of  experimental  configurations.  The  model  con¬ 
stant  evaluation  procedure  is  given  by  Rajendran  and  Grove  [14], 

While  Figure  1  compares  the  model  generated  stress  history  using  the  con¬ 
stants  from  Table  1  with  the  one-dimensional  strain  data  from  the  plate  impact 
test.  Figure  2  compares  the  model  prediction  with  uniaxial  stress  data  from  the 
bar-on-bar  impact  test. 


Modeling  The  Ballistic  Impact  Pressure  Measurements 


Vincent  and  Chang  [18]  conducted  instrumented  ballistic  experiments  on  ceramic 
targets.  A  schematic  of  the  target  configuration  is  shown  in  Figure  3.  Two 
manganin  gauges  were  embedded  into  the  target  assembly:  the  first  (lop)  gauge  is 
placed  between  the  front  ceramic  and  the  isodamp,  and  the  second  (bottom)  gauge  is 
placed  between  the  back  face  of  the  isodamp  and  the  second  ceramic.  This  second 
ceramic  is  backed  by  a  thick  aluminum  plate. 


STEEL 


- -  ■  ■  - 1 -  ■  ■  1 _ _ 

GAUGE#  1  (TOP)v^  AD85  CERAMIC 

AD85  CERAMIC  ^  GAUGE#  2  (BOTTOM) 

ALUMINUM 

Figure  3  A  schematic  of  the  instrumented  ballistic  experiment. 
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The  top  gauge  measurements  seemed  to  he  valid  only  for  a  very  short  time 
(<2  microseconds).  The  ceramic  material  that  surrounds  this  top  gauge  is  destroyed 
by  the  cracked  ceramic.  However,  the  bottom  gauge  survived  during  the  measuring 
period  (about  8  microseconds  to  10  microseconds).  This  gauge  is  protected  by  the 
intact  second  ceramic  and  isodamp.  However,  the  bottom  gauge  is  also  eventually 
destroyed.  Unfortunately,  from  these  destructive  tests  it  is  not  possible  to  determine 
the  timings  of  events  such  as  the  onset  of  microcracking  or  growth  of  macrocracks, 
therefore,  the  different  events  which  might  occur  inside  the  target  during  the  projec¬ 
tile  penetration  can  only  be  indirectly  related  to  the  various  features  of  the  stress 
gauge  signal.  The  data  from  the  two  gauges  are  shown  in  Figure  4.  The  pulse 
duration  and  amplitudes  are  sensitive  to  the  shock  response  of  the  ceramics,  isodamp, 
and  projectile. 


Figure  4.  Manganin  gauge  measured  stress  histories  in  a  ballistic 
experiment 


The  Ballistic  Experiment 

The  ballistic  experiment  of  Vincent  and  Chang  |18|  was  modeled  using  the  86 
version  of  the  EPIC-2  code.  Table  2  provides  the  details  of  simulation  and  experi 
ment.  There  are  six  layers  of  materials  in  the  target  assembly:  (1)  ceramic  (AD85). 
(2)  top  gauge  package,  (3)  isodamp,  (4)  bottom  gauge  package,  (5)  ceramic  (ADC>4). 
and  (6)  aluminum.  The  gauge  records  the  stress-time  history  experienced  by  the 
isodamp.  The  target  dimensions  in  the  simulation  are  given  in  Table  3.  In  the 
simulation,  the  top  and  bottom  ceramics  are  modeled  as  AD85. 


I 


Table  2  Projectile  dimensions  and  materials 

Details 

Experiment 

Modeling 

Length  (mm) 

29  464 

30  0 

Diameter  (mm) 

19  989 

20  0 

Material 

Steel 

Steel 

Velocity 

793  0 

7900 

Table  3.  Target  dimensions 


Layers 

Thickness 

(mm) 

AD85  ceramic 

9  525 

isodamp 

30 

AD85  ceramic 

12  7 

5083  aluminum 

31  75 

Simulation  Details 

The  steel  and  aluminum  were  respectively  modeled  using  the  HY100  steel  and 
2024-T351  aluminum  models  in  the  86  version  of  EPIC-2  library.  Since  the 
aluminum  is  expected  to  remain  mostly  elastic  in  the  present  calculations  due  to 
attentuation  of  the  wave  amplitude,  the  stress  histories  at  the  gauge  locations 
should  not  be  influenced  by  the  plastic  behavior  of  aluminum.  In  fact,  simulations 
with  different  aluminums  produced  almost  identical  stress  histories.  The  Johnson 
and  Cook  strength  model  described  the  strain  rate  dependent  strength  of  these 
materials.  The  Mie-Gruneisen  EOS  was  employed  to  describe  the  bulk  (pressure- 
volume)  behaviors. 

The  isodamp  was  modeled  as  an  elastic-perfectly  plastic  solid  with  a  dynamic- 
yield  strength  of  0.4  GPa  [28].  The  ceramic  strength  model  constants  (see 
Equation  19)  are  Cj  =  4  GPa  and  C3  =  0.03.  The  EOS  is  described  by  the 
Mie-Gruneisen  relationship  (see  Equations  8  and  9).  The  material  density,  shear 
modulus,  and  EOS  constants  for  isodamp  and  AD85  ceramics  are  given  in  Table  4. 
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Table  4.  Material  constants  for  isodamp  and  A085  ceramic 


Material  constants 

A085  ceramic 

Isodamp 

Density  (gm/cm3) 

342 

1  29 

Shear  modulus  (GPa) 

108 

05 

(GPa) 

188 

3.9 

$2  (GPa) 

188 

16  3 

/3a  (GPa) 

0 

420 

r 

1.0 

0  738 

Results  and  Analyses  ^ 

To  understand  the  effects  of  various  deformation  processes  in  the  ceramic  mate¬ 
rial  on  the  calculated  stress  histories  (at  the  gauge  locations),  a  number  of  simula¬ 
tions  of  the  ballistic  experiment  were  performed.  For  this  purpose,  the  following 
cases  were  considered: 

•  Elastic  (Case  E) 

•  Elastic-Plastic  (Case  EP) 

•  Elastic-Cracking  (Case  EC) 

•  Elastic-Plastic-Cracking  (Case  EPC)  » 

These  various  cases  were  simulated  by  properly  adjusting  the  model  parameters. 

For  instance,  the  cracking  is  eliminated  by  setting  the  initial  crack  size  a(,  to  zero. 

The  plastic  flow  is  eliminated  by  setting  the  strength  model  constant  C]  to  a  large 
number.  To  suppress  plastic  pore  collapse,  the  void  content  fu  is  set  to  zero.  In 

the  present  analysis,  the  ramping  of  the  plastic  wave  in  the  plate  impact  test  data  is  *] 

assumed  to  be  due  to  both  strain  rate  and  pore  collapsing  in  AD85.  The  bulk  of 
the  analyses  of  the  EPIC-2  generated  results  is  based  on  the  stress  history  compari¬ 
sons  between  the  experiment  and  each  of  these  cases.  The  effects  of  time  step  and  j 

mesh  on  the  stress  history  were  also  investigated.  i 

I 

» 

Grid  and  Time  Step 

The  EPIC-2  simulation  of  the  ballistic  experiment  idealizes  the  projective  target 
configuration  a:?  an  axisymmetric  geometry,  as  shown  in  Figure  5.  As  shown  in  the 
experiment,  the  radii  of  the  target  layers  in  the  simulations  were  sufficiently  large 
(five  times  the  projectile  radius)  to  prevent  stress  reflections  from  the  lateral  bound- 
aries  influencing  the  calculated  stress  histories.  The  minimum  grid  size  for  the  entire 
mesh  was  selected  based  on  the  number  of  rows  of  elements  and  the  element  size  in 
the  isodamp.  To  ensure  realistic  and  accurate  results,  the  isodamp  was  modeled  using 
four  rows  (layers)  of  elements;  therefore,  for  an  isodamp  layer  of  3  mm  thickness,  the  j 

element  thickness  was  0.75  mm.  This  clement  size  was  the  standard  for  determining  the  p! 

number  of  rows  of  elements  in  the  projectile,  as  well  as  in  the  other  target  layers. 


•  ; 
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The  total  number  of  nodes  in  the  projectile  was  554  and  the  target  was  3888. 
The  number  of  elements  in  the  projectile  was  1040  and  the  target  was  6048.  The 
element  aspect  ratio  was  kept  closer  to  one  in  the  vicinity  of  the  impact  planes  and 
gauge  locations,  as  can  be  seen  from  the  finite  element  grid  shown  in  Figure  6.  In 
general,  the  numerical  results  from  the  shock  wave  based  finite  element/difference 
codes  (called  the  hydrocodes)  arc  inherently  sensitive  to  element  aspect  ratio  and 
the  mesh  type;  therefore,  it  is  important  to  design  suitable  grids  for  realistic  and 
accurate  results  from  the  code  calculations.  For  this  reason,  several  meshes  were 
considered  and  the  repeatability  of  results  was  verified.  While  keeping  the  mesh 
reasonably  fine,  a  variation  in  the  clement  aspect  ratio  was  introduced  between  two 
different  meshes.  The  mesh  with  an  aspect  ratio  of  one  for  the  isodamp  and  two 
for  the  ceramic  compared  very  well  with  the  standard  mesh  with  an  aspect  ratio  of 
one  for  both  isodamp  and  ceramics. 


Figure  6.  A  close-up  view  of  the  finite  element  mesh  near  the  impact  region. 


Figure  7  shows  that  the  calculated  stress  histories  at  the  bottom  gauge  location 
are  compared  for  two  different  meshes.  A  slight  mesh  effect  on  the  stress  histo¬ 
ries  can  be  seen  from  the  plot.  However,  the  overall  shape  and  stress  levels  for 
the  two  meshes  were  practically  the  same.  The  element  size  for  the  ceramic  in 
the  coarse  mesh  was  twice  the  element  size  for  the  isodamp.  Interestingly,  the 
coarse  mesh  with  an  aspect  ratio  of  2  for  the  ceramic  produced  smoother  stress 
history  than  the  fine  mesh. 
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Figure  7.  Effects  of  mesh  on  stress  history 


In  the  EPIC  code,  one  of  the  parameters  that  controls  the  time  step  size  is  the 
”ssf"  parameter.  The  time  step  can  be  controlled  by  setting  values  between  0.1  and 
0.9  for  ssf.  It  is  also  possible  to  examine  the  time  step  effects  through  a  parameter 
which  controls  the  maximum  allowable  lime  step.  The  results  for  the  ssf  =  0.1 
(small  time  step)  and  0.9  (large  time  step)  are  compared  in  Figure  8.  The  stress 
time  histories  are  similar  except  for  some  minor  stress  oscillations.  These  results 
provide  sufficient  confidence  in  the  numerical  results.  One  of  several  other  code 
parameters,  such  as  the  maximum  allowable  time  step,  can  also  influence  the  numeri¬ 
cal  results;  however,  the  proper  choice  of  these  time  step  related  code  parameters 
produce  similar  and  repeatable  results. 

Elastic  (Case  E) 

The  simplest  stress-strain  relationship  is  the  Hooke’s  Law  for  an  elastic  material. 
Since  ceramic  is  a  brittle  solid,  it  is  proper  to  begin  the  analysis  with  an  elastic 
description.  The  inelastic  strains  due  to  both  microcracking  and  plastic  flows  are 
not  allowed  in  the  elastic  simulation.  The  shear  and  bulk  moduli  will  not  degrade 
and  will  remain  intact  under  impact  loading.  The  strength  of  the  ceramic  is  unlim¬ 
ited;  therefore,  failure  is  not  allowed  in  the  elastic  case.  The  results  from  this  case 
is  presented  in  Figures  9  and  10.  The  bottom  stress  gauge  data  arc  compared  with 
the  simulation  in  Figure  9.  As  can  be  seen  from  this  figure,  the  calculated  stress 
levels  are  higher  and  the  loading  duration  is  lower  when  compared  to  the  data. 

These  results  clearly  indicate  that  a  simple  elastic  assumption  is  inadequate  for 
describing  the  complex  impact  behavior  of  ceramic  materials. 
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Figure  9.  A  comparison  between  bottom  gauge  data  and  Case  E 
The  ceramic  behavior  is  assumed  elastic. 
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Figure  10.  A  comparison  between  top  gauge  data  and  Case  E  The 
ceramic  behavior  is  assumed  elastic 


A  comparison  between  the  top  gauge  and  the  elastic  simulation,  as  shown 
in  Figure  10,  further  confirms  the  inadequacy  of  simple  elastic  description  of  the 
ceramic.  However,  the  top  gauge  record  is  questionable  beyond  2  microseconds. 

In  the  experiment,  the  fractured  ceramic  that  surrounds  the  gauge  destroyed  the 
gauge  beyond  this  time.  This  was  confirmed  by  Vincent  and  Chang  |1X|  through 
several  experiments.  The  gauge  measurements  did  not  show  the  arrival  of  a  second 
strong  shock  as  it  did  in  the  simulation;  the  simulation  clearly  showed  the  arrival  of  a 
strong  second  shock.  A  relatively  simple  wave  analysis  under  one-dimensional  strain 
condition  (as  in  the  plate  impact  experiment)  shows  that  the  second  shock  arrives 
from  the  boundary  between  the  isodamp  and  the  bottom  ceramic  due  to  an  over¬ 
whelming  impedance  mismatch  between  the  two  layers. 

The  absence  of  such  second  shock  in  the  test  imply  that  the  top  gauge  data  ma\ 
be  reliable  only  for  the  first  2  microseconds.  A  carbon  gauge  (places  in  addition  to 
a  manganin  gauge)  showed  the  arrival  of  a  second  shock;  however,  this  gauge  was 
also  destroyed  within  3  microseconds  after  the  impact.  It  is  also  possible  that  the 
apparent  second  shock  in  the  carbon  gauge  signal  could  be  an  artifact  of  gauge  fail¬ 
ure.  Since  only  one  experiment  was  conducted  with  the  carbon  gauge,  a  definite  con¬ 
clusion  could  not  be  made.  The  measured  maximum  stress  amplitude  in  the  top 
gauge  seems  to  be  a  reliable  data.  An  approximate  one -dimensional  strain  analysis 
based  on  the  impact  shock  amplitude  in  the  ceramic  and  the  impedance  match  solu¬ 
tion  between  the  isodamp  and  ceramic  seems  to  support  this  conclusion. 
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Elastic-Plastic  (Case  EP) 


In  the  elastic-plastic  simulation,  the  impact  behavior  of  ceramic  was  described 
through  the  strain  rate  dependent  strength  in  Equation  1^.  An  initial  value  ol  zero  • 

for  the  microcrack  size  eliminates  microcracking  in  the  calculation  The  initial  poros¬ 
ity  was  also  set  to  zero  so  that  pore  collapse  would  not  occur.  Since  the  shock 
waves  in  the  ballistic  experiment  are  spherical,  it  is  not  possible  to  accurately  deter 

mine  the  different  wave  arrival  times.  However,  the  rough  estimates  of  shock  arrival  j 

times  based  on  the  one-dimensional  strain  analysis  sometimes  provide  guidelines  lor  j 

interpreting  the  gauge  signal.  In  Figure  11,  the  elastic  wave  arrives  at  the  gauge  loca¬ 
tion  at  about  2  microseconds.  A  weaker  release  wave  from  the  edge  of  the  projec¬ 
tile  follows  this  initial  elastic  wave  and  unloads  at  point  A.  A  reshock  (compressive 
loading)  later  arrives  from  the  interface  of  the  isodamp  and  bottom  ceramic  at  point 
B.  The  unloading  waves  from  the  free  surfaces  of  the  top  ceramic  unload  the 

stresses  at  point  C.  In  the  experimental  data,  the  sharp  peak  at  point  A.  as  well  as  •  : 

the  strong  second  shock,  are  not  found. 
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Figure  1 1  A  comparison  between  bottom  gauge  data  and  Case  EP 
The  ceramic  behavior  is  assumed  elastic-plastic  (no  cracking) 


When  the  ceramic  behavior  is  described  by  an  elastic-plastic  model,  the  model 
comparison  with  the  bottom  gauge  data  improved  significantly,  as  shown  in  Figure  II 
It  appears  that  limiting  the  ceramic  compressive  strength  to  finite  values  through  a 
yield  surface  could  improve  the  model  prediction.  With  lack  of  anv  microstructural 
evidence  to  prove  macroplaslic  flows  in  the  brittle  ceramic  under  impact  loading  condi 
tions,  it  is  premature  to  conclude  that  ADN5  ceramic  deforms  plastically  like  a  metal 


» 
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just  because  the  clastic-plastic  model  reproduced  the  experimental  measurements. 

There  is  also  a  possibility  to  cap,  or  limit,  the  ceramic  strength  due  to  microcrackmg 
and  crushing.  Though  the  matching  is  good,  the  absence  of  certain  salient  features 
in  the  simulation  indicates  elastic-plastic  idealization  alone  may  not  be  reproducing  all 
the  features.  It  appears  that  the  matching  between  the  simulation  and  experiment 
can  be  significantly  improved  by  limiting  the  ceramic  strength  to  finite  values  (see 
Equation  19).  To  further  verify  the  results  from  the  elastic-plastic  idealization,  a  compar¬ 
ison  between  simulation  and  top  stress  gauge  data  is  made  in  Figure  12.  The  arrival  cl' 
the  second  shock  at  point  B  can  be  clearly  seen  in  the  simulated  stress  history'. 


Figure  12  A  comparison  between  top  gauge  data  and  Case  EP  The  ceramic 
behavior  is  assumed  elastic-plastic  (no  cracking) 

Elastic-Cracking  (Case  EC) 

This  case  examines  the  effect  of  microcracking  on  the  stress  profile.  In  the 
simulation,  a  large  value  for  the  constant  C|  in  Equation  19  was  assumed  to  elimi¬ 
nate  plastic  flows.  The  strength  was  degraded  or  relaxed  through  microcracking  only. 
The  viscoplastic  pore  collapsing  was  also  suppressed  in  the  simulation.  Figure  In  com¬ 
pares  the  calculated  stress  history  from  the  elastic-cracking  case  and  the  bottom  gauge 
data.  The  two  stress  peaks  (at  points  A  and  C)  wcic  present  as  irt  Cases  E  and  EP. 


The  stress  levels  are  instantaneously  higher  compared  to  the  gauge  data.  The  stress 
history  during  loading  and  unloading  matched  with  the  data  somewhat  in  an  average 
manner.  Interestingly,  the  results  from  Case  EP  (elastic-plastic)  compared  relatively 

better  with  the  data  than  this  elastic-cracking  case.  These  results  further  confirm  § 

that  the  ceramic  behavior  is  much  more  complex  than  the  assumptions  of  either 
elastic-cracking  or  elastic-plastic. 
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Figure  13  A  comparison  between  bottom  gauge  data  and  Case  EC  The  ceramic 
behavior  is  assumed  elastic  cracking  (no  plastic  flow) 

» 

The  Full  Ceramic  Model  (Case  EPC) 

To  further  investigate  the  effects  of  inelastic  deformations  on  the  stress  history, 
the  elastic-plastic  cracking  case  with  pore  collapse  was  considered.  In  other  words.  • 

the  Rajendran-Grove  ceramic  model  was  used  to  describe  the  impact  behavior  of 
AD85  under  Case  EPC.  The  ballistic  impact  experiment  of  Vincent  and  Chang  |1<S| 
was  simulated  using  the  ceramic  model  constants  in  Table  1.  A  10%  porosity  content 
and  an  initial  Haw  size  of  0.057  mm  were  used  in  the  simulation.  Figure  14 
compares  the  model  and  the  bottom  gauge  data. 
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Figure  14.  A  comparison  between  the  ceramic  model  (with  tensile 
ni  =  1)  and  bottom  gauge  data. 


The  model  showed  a  lower  stress  level  and  a  higher  pulse  duration  compared  to 
the  data.  The  analysis  indicated  an  excessive  tensile  damage  in  the  ceramic.  Several 
regions  of  the  top  ceramic  plate  had  fractured.  Since  the  Rajendran -Grove 
ceramic  model  degrades  the  ceramic  strength  due  to  microcracking.  the  microcracks 
induced  damage  significantly  lowering  the  stress  amplitude. 

To  further  investigate  the  ceramic  model  results,  damage  shade  plots,  as  shown  in 
Figures  15  and  16,  were  generated.  The  regions  without  any  shade  are  free  of  dam¬ 
age.  The  damage  evolution  at  4  microseconds  is  shown  in  Figure  15.  The  maximum 
damage  was  about  0.14.  This  occurred  at  a  location  2  mm  in  depth  and  a  13  mm 
radius  (from  the  axis  of  symmetry)  which  corresponds  to  the  Hertzian  crack  vicinity. 
The  ceramic  model  could  indicate  the  formation  of  classical  Hertzian  cracks  emanat¬ 
ing  at  about  a  45°  angle  from  the  edge  of  the  projectile-target  interaction  regions. 

The  damage  shade  plot  at  8  microseconds  is  shown  in  Figure  l(>.  The  maximum 
damage  occurred  in  regions  closer  to  the  top  gauge  location.  Since  the  isodamp  is  a 
low  impedance  and  low  strength  material,  the  impact  loading  conditions  create 
large  biaxial  tensile  stresses  in  the  ceramic.  Therefore,  at  the  gauge  location,  tensile 
damage  evolves  rapidly  and  fractures  the  ceramic.  Since  ceramic  is  extremely  weak  in 
tension  and  strong  under  compressive  loading,  the  tensile  regions  start  fracturing  first. 
When  damage  reaches  a  value  of  one.  the  ceramic  is  assumed  to  have  failed  completely, 
and  the  failed  elements  will  not  sustain  any  tensile  loading.  However,  those  elements 
continue  to  carry  compressive  loading  until  the  material  pulverizes  under  compression. 


Figure  16  The  microcracking  damage  shade  ptot  at  time  8  microseconds. 


While  the  projectile  penetrates  into  the  top  ceramic  layer,  damage  due  to 
microcracking  develops  in  the  bottom  ceramic  which  is  separated  from  the  lop 
ceramic  layer  by  the  isodamp.  The  plots  of  damage  shade  in  the  bottom  ceramic 
layer  corresponds  to  very  low  values  (<0.05)  of  the  damage.  The  mild  plastic  How 
near  the  axis-of-symmetry  at  the  top  surface  of  the  aluminum  causes  tensile  loading 
in  ceramic.  These  tensile  stresses  initiate  damage  in  the  bottom  ceramic,  as  shown  in 
Figure  16. 

In  the  analysis,  a  value  of  1.0  for  n|  overestimated  the  tensile  damage  in  the  top 
ceramic  plate.  According  to  the  microcrack  growth  rate  in  Equation  16,  the  limiting 
crack  growth  rate  is  equal  to  the  Rayleigh  wave  speed  when  n^  =  1;  therefore,  to 
improve  the  matching  between  the  experiment  and  the  model,  the  tensile  crack 
growth  index  n|  is  reduced  from  1.0  to  0.1  and  the  ballistic  experiment  was  simulated 
again.  All  the  other  constants  in  Table  1  were  not  modified.  Recall  that  both  the 
elastic  and  elastic-plastic  simulations  showed  a  strong  second  shock  arriving  at  the  top 
gauge  location  at  about  2.4  microseconds.  Though  the  data  did  not  show  any  second 
shock,  the  maximum  amplitude  of  the  first  shock  recorded  by  the  gauge  was  a  reli¬ 
able  data.  Therefore,  to  further  verify  the  stress  history  at  the  top  gauge  location,  a 
comparison  between  the  model  and  the  data  is  made  in  Figure  17.  As  can  be 
seen  from  this  figure,  the  elastic-plastic  cracking  simulation  also  showed  the  second 
shock  at  about  2.4  microseconds. 
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Figure  17  A  comparison  between  top  gauge  data  and  the  ceramic  model 
(with  tensile  m  =  0.1). 


It  is  also  worth  recalling  that  the  simulated  peaks  in  the  elastic  and  elaslic-plaslie 
cases  were  40  Kbars  and  30  Kbars  (10  Kbars  =  1  GPa),  respectively.  These  values 
are  significantly  higher  compared  to  the  measured  maximum  stress  of  19  Kbars.  How¬ 
ever,  the  model  calculated  peak  of  about  17  Kbars  compares  well  with  the  data. 

To  further  analyze  the  ceramic  model  predictions,  a  comparison  between  the 
model  generated  stress  history  and  measured  stress  history  (bottom  gauge)  is  shown  in 
Figure  18.  As  can  be  seen  from  this  figure,  the  model  reproduced  the  experimental 
data  extremely  well.  Both  the  amplitude  and  the  pulse  duration  matched  between  the 
model  and  the  data.  The  value  of  0.1  for  the  tensile  crack  growth  rate  parameter 
improved  the  model  prediction  significantly.  Since  this  parameter  indirectly  controls 
the  amount  of  stress  relaxation  in  the  model,  a  reduced  value  for  this  parameter  ac¬ 
cordingly  increased  the  calculated  stress  amplitude. 


Figure  18.  A  comparison  between  the  Rajendran-Grove  ceramic 
model  (with  tensile  n»  =  0. 1)  and  bottom  gauge  data. 


The  results  from  these  simulations  indicate  that  the  model  prediction  based  on 
both  brittle  microcracking  and  plastic  How  in  the  ceramic  material  matched  the  experi¬ 
mental  measurements  extremely  well.  The  ceramic  model  not  only  matched  the  stress 
amplitude  and  the  time  duration,  it  also  reproduced  most  of  the  salient  features  of 
the  measured  stress  signal.  The  shade  plots  of  damage  also  showed  the  experimen¬ 
tally  observed  fracture  pattern  such  as  the  Hertzian  cracks  emanating  from  the  edges  of 
projectile-target  interaction  regions  and  the  ceramic  fracturing  around  the  top  gauge. 
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As  a  final  exercise,  the  effect  of  interfaces  between  the  target  layers  (top  ceramic 
isodamp,  isodamp  bottom  ceramic,  and  bottom  ceramic  aluminum)  on  the  calculated 
stress  history  was  investigated.  Two  cases  were  considered:  (1)  with  slide  lines 
between  the  layers,  and  (2)  without  any  slide  lines.  In  case  (1)  the  interlaces  are 
allowed  to  slip,  and  in  case  (2)  interfaces  are  glued  together  so  that  no  slip  is 
allowed.  The  results  from  these  two  cases  are  compared  in  Figure  19. 


Figure  19.  A  comparison  between  stress  histories  for  the  cases  with  and  without  interface 
slide  lines  and  the  bottom  gauge  data. 


The  comparison  plot  in  Figure  19  clearly  shows  that  the  calculated  stress  histories 
with  and  without  slide  lines  produced  similar  results.  This  comparison  demonstrates 
that  the  effect  of  interface  sliding  on  the  stress  history  is  not  significant,  at  least  in 
the  present  application.  However,  the  slide  lines  introduced  some  sort  of  (periodic) 
oscillations  in  the  stress  history  during  unloading.  This  may  be  due  to  the  slide-line 
interaction  introduced  numerical  noise.  These  results  show  that  the  influence  of  time 
step,  mesh  size,  and  slide  line  on  the  model  calculated  stress  histories  is  minimum. 
Therefore,  it  may  be  concluded  that  the  differences  between  the  stress  histories  for 
the  various  cases  are  truly  due  to  different  types  of  material  behaviors. 
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Summary  and  Conclusions 


Summary 

Recently,  Rajendran  [11]  reported  the  development  of  an  advanced  ceramic 
model.  This  Rajendran-Grove  ceramic  model  has  been  implemented  into  the  86  ver¬ 
sion  of  the  EPIC-2  code  and  successfully  used  to  model  the  impact  behavior  of  AD85 
ceramic  under  different  impact  test  configurations  such  as  the  plane  plate  impact  and 
rod-on-rod  tests.  The  stress-strain  states  under  these  two  configurations  are  fairly  sim¬ 
ple.  In  the  plastic  impact  test,  the  strain  is  one-dimensional,  and  in  the  rod-on-rod 
test  the  stress  is  one-dimensional.  There  are  seven  model  constants: 

Kjc,^,  N*.  ao,  n  j  ,  n  J"  ,  and  n2-  The  fracture  toughness  value  is  taken  from  the 
fracture  mechanics  handbooks.  The  rest  of  the  constants  were  determined  from  the 
experimental  data.  For  this  purpose,  the  manganin  gauge  measured  stress  histories 
from  the  plate  impact  and  rod-on-rod  impact  experiments  were  employed. 

The  main  objective  of  the  present  work  is  to  demonstrate  the  applicability  and 
generality  of  the  Rajendran-Grove  ceramic  model  under  relatively  complex  stress-strain 
states.  The  idea  is  to  employ  the  AD85  constants,  determined  from  the  plate  impact 
and  rod-on-rod  impact  tests  data,  to  successfully  describe  the  AD85  ceramic  behavior 
under  an  entirely  different  experimental  configuration.  This  report  presented  the 
ceramic  model  in  detail  and  its  successful  application  to  describe  the  deformation  and 
fracture  in  a  layered  ceramic  target  due  to  a  steel  projectile  impact.  The  three- 
dimensional  stress-strain  state  under  this  target  configuration  is  fairly  complex  due  to 
shock  wave  interactions.  In  the  ballistic  experiment,  two  manganin  gauges  were 
embedded  in  the  target.  This  ballistic  experiment  was  conducted  by  Vincent  and 
Chang  [18]. 

The  86  version  of  the  EPIC-2  code  was  used  to  simulate  the  ballistic  experiment. 
Several  simulations  were  performed  to  understand  and  evaluate  the  effects  of  elastic, 
elastic-plastic,  elastic-cracking,  and  elastic-plastic  cracking  behaviors  on  the  measured 
stress  histories.  The  results  from  the  simulation  in  which  the  ceramic  is  assumed  to 
behave  fully  elastically  (without  plastic  flow  and  microcracking)  did  not  match  the 
data,  whereas  the  elastic-plastic  assumption  produced  reasonable  results.  The  overall 
matching  between  the  data  and  simulation  was  reasonably  good.  However,  the  details 
of  the  simulated  wave  profile  did  not  agree  with  the  experimentally  measured  stress 
profile.  Similar  results  were  observed  when  the  ceramic  was  modeled  as  an  elastic¬ 
microcracking  solid  without  any  plastic  flow.  These  results  showed  that  either  the 
elastic-plastic  or  the  elastic-cracking  behavior  was  adequate  to  describe  the  impact 
behavior  of  ceramic.  To  verify  the  effects  of  a  combined  elastic-plastic  cracking 
behavior  on  the  stress  history,  the  Rajendran-Grove  ceramic  model  was  employed  in 
the  ballistic  experimental  simulation. 

In  the  simulation,  the  AD85  constants  in  Table  1  (with  the  exception  of  the  ten¬ 
sile  crack  growth  index  constant,  nj)  were  employed  to  describe  the  ceramic  as  an 
elastic-plastic  microcracking  material.  The  calculated  stress  history  matched  the  mea¬ 
sured  stress  history  data  extremely  well.  When  a  value  of  1  was  employed  for  n^. 
the  ceramic  exhibited  extensive  damage  and  the  simulated  stress  profile  did  not  match 
with  the  experiment:  however,  a  value  of  0.1  for  nj  predicted  the  data  successfully. 
The  Rajendran-Grove  model  also  showed  physically  possible  fracture  patterns  in  the 
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ceramics.  In  summary,  the  ceramic  model  constants,  estimated  from  the  standard 
one-dimensional  impact  tests  data,  reproduced  the  multiaxial  ballistic  experimental  data 
well. 

Conclusions 

The  inelastic  deformations  in  ceramic  materials  due  to  impact  loading,  in  general, 
consists  of  elastic,  plastic,  and  microcracking  components.  It  is  not  experimentally 
possible  to  isolate  these  deformations  through  any  direct  measurements.  With  lack 
of  any  recovery  techniques  to  examine  the  post-impacted  ceramic  targets,  especially  at 
high  velocity  impact,  only  speculative  assumptions  can  be  made  on  the  various  defor¬ 
mation  and  fracture  processes. 

The  computational  analysis  of  the  impact  experiments  using  advanced  ceramic 
models  will  indeed  help  in  evaluating  the  various  possible  deformation  and  fracture 
modes  in  the  ceramic  materials.  So  far,  the  experimentalist  have  been  making  accu 
rate  velocity  and  stress  measurements,  high  speed  photographs,  and  X-ray  radiographs. 
Interpretations  and  validations  of  these  valuable  measurements  demand  a  detailed  com 
putational/analytical  modeling  of  the  impact  experiments.  These  modeling  efforts  will 
eventually  lead  to  a  greater  understanding  of  the  impact  behavior  of  ceramics. 

In  order  to  increase  confidence  in  using  the  ceramic  models  as  predictive  tools 
in  armor/antiarmor  applications,  generality  of  the  model  parameters  should  be  tested 
using  a  variety  of  experimental  data.  The  experiments  should  not  only  include  vuri 
ous  stress-strain  states,  but  also  a  range  of  velocity  regimes.  Though  the  ballistic 
experiment  could  model  well  using  the  constants  determined  from  one-dimensional 
experiments,  the  stress-strain  state  dependency  on  the  tensile  cracking  requires  addi¬ 
tional  investigation;  therefore,  the  effects  of  the  initial  microflaw  size  and  the  number 
of  microflaws  on  microcracking  should  be  studied  through  additional  simulations.  In 
general,  the  continuum  mechanics  based,  three-dimensional  Rajendran-Grove 
ceramic  model  produced  realistic  results  and  matched  a  wide  variety  of  experimental 
impact  data. 
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Appendix 


Degraded  Moduli  Expressions 

Margolin  [15]  derived  the  degraded  elastic  moduli  for  noninteracting,  penny-shaped 
microcracks  of  various  sizes  and  in  random  orientations  and  the  corresponding  expres¬ 
sions  for  the  isotropic  elastic  moduli,  in  Equation  3  are: 

Qjkl  =  Cj.dfcdji  +  C2 <5jk <5j|  +  C3<5ik<5j|  ,  (19) 

where 

C]  =  Bq  +  ,  (20) 


Cs  =  ^  2(1  +  v)G  ■  (22) 

In  the  above  equations,  G  and  v  are  the  shear  modulus  and  Poisson's  ratio,  respec¬ 
tively,  of  the  undamaged  material  while  A0,  B0,  and  D()  are  damage  parameters  whose 
values  depend  on  the  stress  state.  To  evaluate  these  parameters,  Margolin  defined  a 
microcrack  density  parameter. 


a 


3 

max 


45  E 


(23) 


In  Equation  23,  y0*  is  the  microcrack  density  and  E  is  the  Young's  modulus  of  the  un¬ 
damaged  material.  N*  is  the  number  of  microcracks  per  unit  volume,  and  amax  is  the 
maximum  microcrack  size.  Margolin  identified  the  following  four  cases  of  stress  state 
in  evaluating  the  damage  parameters  A(>  Bt).  and  D0: 

Case  1:  o\ ,  <72 ,  <73  >  0 

(all  principal  stresses  are  tensile) 
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Ao  =  (1  -  vz)  -  (1  +  »)]  y 


(24) 


*o  =  [(1  -  v2)  +  4(1  +  v)  y 


D0  =  Ao 


Case  2:  a\ ,  ,  a->,  <  0 

(all  principal  stresses  are  compressive) 


Ao  =  -  (1  +  y)y* 


B0  =  4(1  +  v)y* 


D0  —  Aq 


Equations  27  through  29  will  result  in  the  degradation  of  the  shear  modulus,  but  not 
of  the  bulk  modulus  because  under  compression  only  crack  movement  of  the  closed 
microcracks  under  modes  II  and  III  are  permitted. 


Case  3:  at .  <72 ,  >  0 .  <73  <  0 

(two  principal  stresses  are  tensile  and  one  principal  stress  is  compressive) 


Ao  — 


93  _  x 


(1  -  O  ~  (1  +v)\  y* 


9-3  _  -j/jSx 


2  ^  >  (1  -  vz)  +  4(1  +  v)  y* 


D0  =  [(6^5  -  5/J3)  (1  -v2)  -  (!+«;)]  7' 


where 
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fi 


°\  +  a2 


o  |  +  02  ~  2  <73 


(33) 


Case  4:  ctj  >  0 ,  <72 ,  <73  <  0 

(one  principal  stress  is  tensile  and  two  principal  stresses  are  compressive) 


Ao  — 


Bo  = 


5(1  -  3(1  -ft  (]  ,.2).  (1+y) 


^-A-3^5)  (1.w2)  +  4(1+|;) 


D0  =  {[6(1  -/?5)  -  5(1  -y53)]  (1  -v2)  -  (1  +1.)}  y 


(34) 


(35) 


(36) 


where 


P 


<72  +  <73 


<72  +  <73  —  2  <7 1 


(37) 


In  the  Rajendran-Grove  ceramic  model,  the  compliance  tensor  C  is  analytically  in¬ 
verted  to  the  stiffness  tensor  M  using  the  following  identity  relationship: 


_  (^ik  <*il  +  <*il  ^jk) 

Mijkl  Cjjki - 2 - ; 


(38) 


a  relationship  can  be  established  between  C\  C2.  and  C3  of  Equation  20  through  22 
and  mi,  m2,  and  m3  of  the  following  relationship: 


M,jki  =  m!  <3lk  <3j!  +  m2  djjrSjk  4  m3  <3,j<5k| 


(39) 


where 


m  1  =  m2 


2(C[  +  C  2) 


(40) 


and 
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m3  - 


•Cy _ _ 

-s  +  c2  +  V',) 


The  effective  degraded  shear  and  bulk  moduli  are  defined  as 


G  —  mi 


m3  +  j mi 


Instead  of  employing  the  Margolin's  expression  for  Case  1  (all  principal  stresses  are 
tensile),  the  Budiansky  and  O’Connell  [20]  solutions  can  be  used  for  randomly  ori¬ 
ented,  noninteracting  microcracks  under  tensile  loading.  The  corresponding  relation¬ 
ships  for  the  damaged  stiffness  solutions  are: 


mj  =  m2  =  G 
2  G  v 

m3  =  - — 

(t-2v  ) 


where 


1  16 

v  =  v  1  -  33-y 


G  =  G  1 


32)  0  ~  i'X5  “  ^ 

45  j  (2  -  v)  7 


In  these  equations,  v  and  G  are  the  Poisson's  ratio  and  shear  modulus,  respectively,  ol 
the  microcrack  damaged  material.  Using  m|,  and  m3  from  Equations  44  and  45. 
the  degraded  bulk  modulus  can  be  computed  from  Equation  43.  It  is  obvious  from 
Equations  44  through  47  that  a  complete  loss  of  strength  is  predicted  when  the 
microcrack  density  y  reaches  9/16.  For  tensile  loading  conditions,  based  on  the  com¬ 
parison  with  Margolin's  equations,  there  is  no  bound  on  the  crack  density.  However, 
Budiansky  and  O’Connell’s  solution  limits  the  crack  density  to  9/16.  This  permits  the 
damage  parameter  to  vary  from  zero  (no  damage)  to  one  (fully  damaged).  Therefore,  in 
the  ceramic  model,  Budiansky  and  O’Connell’s  equations  are  used  instead  of  Margolin's 
equations  for  the  case  when  all  the  principal  stresses  are  positive  (see  Case  1). 
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